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The geometric frustration in a class of the mixed spin-1/2 and spin-S Ising-Heisenberg diamond chains is 
investigated by combining three exact analytical techniques: Kambe projection method, decoration-iteration 
transformation and transfer-matrix method. The ground state, the magnetization process and the specific 
heat as a function of the external magnetic field are particularly examined for different strengths of the 
geometric frustration. It is shown that the increase of the Heisenberg spin value S raises the number of 
intermediate magnetization plateaux, w/hich emerge in magnetization curves provided that the ground state is 
highly degenerate on behalf of a sufficiently strong geometric frustration. On the other hand, all intermediate 
magnetization plateaux merge into a linear magnetization versus magnetic field dependence in the limit of 
classical Heisenberg spin S —> oo. The enhanced magnetocaloric effect with cooling rate exceeding the one 
of paramagnetic salts is also detected when the disordered frustrated phase constitutes the ground state and 
the external magnetic field is small enough. 
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1. Introduction 

The quantum Heisenberg model with diamond chain topology has enjoyed a great scientific 
interest since two unusual tetramer-dimer and dimer-monomer phases were theoretically predicted 
by Takano et al. [ [1] in the zero-field ground-state phase diagram of the spin-1/2 Heisenberg 
diamond chain as a result of the mutual interplay between quantum fluctuations and geometric 
frustration. Motivated by this discovery, several other one-dimensional (ID) quantum spin models 
consisting of diamond-shaped units have been suggested and solved with the aim to bring insight 
into a frustrated magnetism of diamond chain systems. It is now widely accepted that the zero- 
field ground-state phase diagram of the spin-1/2 Heisenberg model with the distorted diamond 
chain topology is rather complex and consists of the usual ferrimagnetic phase as well as several 
quantum dimerized and plaquette states [ [H [31 [H [SI [5] . Further theoretical studies devoted to 
this ID quantum spin model have provided accurate results for the ground-state phase diagram 
in a presence of the external magnetic field [ 13 [H [71 [8] , the spin gap [ [9] , the magnetization 
and susceptibility [[10], as well as, the inversion phenomenon, which can be induced through the 
exchange anisotropy [ [Tl] [12] . It is noteworthy that the ground state and thermodynamics of the 
mixed-spin diamond chains constituted by higher spins have been particularly examined as well [ 
[ni[ll[15l[16l[17]. 

It should be also mentioned, however, that the immense theoretical interest focused on the 
diamond chain structures is not purposeless. An important stimulus for a theoretical treatment 
of diamond chain models comes from rather recent findings that several insulating magnetic ma- 
terials such as azurite Cu3(C03)2(OH)2 [ [H [ISl [201 [211 [22l [23l [24] , Bi4Cu3V20i4 [ l25j, and 
Cu3(Te03)2Br2 [H^, represent possible experimental reahzations of the spin-1/2 Heisenberg di- 
amond chain. It is convenient to remark that the natural mineral azurite was for a long time 
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regarded as the best known candidate of the diamond chain compound even though recent ab- 
initio calculations indicate possibly non-negligible intcr-chain interactions in this compound [i27j. 
Notwithstanding this fact, experimental data measured for the azurite are in a good qualitative 
accordance with the relevant theoretical predictions for the highly frustrated spin- 1/2 Heisenberg 
diamond chain. As a matter of fact, the magnetization plateau at one-third of the saturation 
magnetization [ [191 EOl EH [22l [23] , the double-peak structure in temperature dependences of the 
specific heat [ HSl |H [H [H] and the zero-field susceptibility [ [H [TH HOI HI [13] have been 
found both theoretically as well as experimentally. In addition, the spin-1, spin-3/2, and spin- 
5/2 Heisenberg model with the diamond chain topology might prove its usefulness in elucidating 
magnetic properties of polymeric coordination compounds M3(OH)2 (M = Ni, Co, Mn) [ \28\ I29j. 
[Ni3(fum)2(M3-OH)2(H20)4]-(2H20) [^J and [Co3(NC5H3(COO)2)2(Ai3-H20)2(H20)2] [[31]. 

Unfortunately, the rigorous theoretical treatment of geometrically frustrated quantum Heisen- 
berg models is very difficult to deal with due to a non-commutability of spin operators involved in 
the Heisenberg Hamiltonian, which is also a primary cause of a presence of quantum fluctuations. 
Owing to this fact, we have recently proposed a novel class of the geometrically frustrated Ising- 
Heisenberg diamond chain models [ |32l [33l [Ml [35] , which overcome this mathematical difficulty 
by introducing the Ising spins at the nodal sites and the Heisenberg dimers on the interstitial 
decorating sites of the diamond chain. This simplified quantum model can be examined within the 
framework of exact analytical approach based on the generalized decoration-iteration transforma- 
tion [ 36, i37l 138] , because the nodal Ising spins represent a barrier for quantum fluctuations that 
are consequently restricted to elementary diamond-shaped units. It is worth mentioning that the 
relatively simple analytical technique based on the generalized decoration-iteration transformation 
has been recently adapted to explore an effect of the geometric frustration also in the asymmetric 
Ising-Heisenberg tetrahedral chain [[3S], the Ising- Heisenberg chain consisting of triangular-shaped 
Heisenberg trimers alternating with the nodal Ising spins [|40l[4T], as well as, the kinetically frus- 
trated diamond chain models constituted by the nodal Ising spins and mobile electrons delocalized 
over the interstitial decorating sites [ |42l [43l [44] . 

The main purpose of this work is to provide the exact solution for the generalized version of 
the mixed spin-1/2 and spin-S* Ising-Heisenberg diamond chain [ [3H [331 [32 [3S], which should 
bring a deeper insight into how the magnetic properties depend on the quantum spin number 
S of the Heisenberg spins. The exact analytical solution for this extended version of the Ising- 
Heisenberg diamond chain will be attained by combining the Kanibe projection method [[45] with 
the generalized decoration- iteration mapping transformation [ 1361 137|, 138] and the transfer-matrix 
technique [[iSlITT]. 

The outline of this paper is as follows. In section [2l we present at first a detailed formulation 
of the Ising-Heisenberg model, which is subsequently followed by a brief description of basic steps 
of exact analytical treatment. Section [3] deals with the interpretation of the most interesting re- 
sults for the ground-state phase diagrams, the magnetization process, the specific heat, and the 
magnetocaloric effect. Finally, some concluding remarks are drawn in section [2 

2. Model and its exact solution 

Let us consider an one-dimensional lattice of inter-connected diamonds as schematically illus- 
trated in figure m In this fi gure, empty circles denote nodal lattice sites occupied by the Ising 
spins /i — 1/2, while the filled ones label interstitial (decorating) lattice sites occupied by the 
Heisenberg spins with an arbitrary quantum spin number S. The total Hamiltonian of this mixed- 
spin Ising-Heisenberg diamond chain can be for further convenience written as a sum over cluster 
Hamiltonians 7i = X^a^i '^k, where each cluster Hamiltonian Tifc involves all the interaction terms 
belonging to the fcth diamond-shaped unit (see figure [1]) 

Hk = JH^3k-r^-ik+Ji{Slk-i+SI,k){i^lk-2+i^lk+i)-HH{Slk-i+S^k)^Hi{ff^^^ (1) 

Here, /i,^ and — {S'^, S^, S^) denote spatial components of the spin-1/2 and spin-S" operators, the 
parameter Jh labels the isotropic exchange interaction between the nearest-neighbouring Heisen- 
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berg spins and the parameter J/ denotes the Ising interaction between the Heisenberg spins and 
their nearest Ising neighbours. FinaUy, the last two terms determine the magnetostatic Zeeman's 
energy of the Ising and Heisenberg spins placed in an external magnetic field Hj and Hh oriented 
along the z-axis, respectively. 




^h unit 

Figure 1. A part of the mixed spin-1/2 and spin-S Ising-Heisenberg diamond chain. The empty 
circles denote lattice positions of the Ising spins /i = 1/2, while the filled circles label lattice 
positions of the Heisenberg spins of an arbitrary magnitude 5*. The ellipse demarcates spins 
belonging to the fcth diamond unit. 

The most important point of our calculation represents an evaluation of the partition function. 
Taking into account a validity of the commutation relation \Hk,'Hi\ = between cluster Hamilto- 
nians of two different diamond units, the partition function of the Ising-Heisenberg model can be 
partially factorized into the product 

N 

Z^Y. Tr{5.} exp(-/37T;) =Y.X{ Trfc exp(-/3Hfc)- (2) 

Above, (3 = l/(fcBT), fee is Boltzmann's constant, T is the absolute temperature, the symbols 
and Tr{5^} denote summations over spin degrees of freedom of all Ising and Heisenberg 
spins, respectively, and the symbol Tr^ stands for a trace over spin degrees of freedom of both 
Heisenberg spins from the kth diamond plaquette. It is quite obvious from the equation ^ that 
it is necessary to perform this latter partial trace in order to proceed further with a calculation. 
The calculation of the partial trace can be easily accomplished with the help of Kambe projection 
method [ 45 , since the cluster Hamiltonian Ti.^ can alternatively be viewed as the Hamiltonian 
of the spin-S* Heisenberg dimer placed in the effective field H^ff = Hh — Ji{f^lk-2 + /^3/c+i)- 
Consequently, the complete set of eigenvalues Ek corresponding to the cluster Hamiltonian ([1]) can 
be expressed solely in terms of two quantum spin numbers Stot and S^qj, 

EkiStot, S^t) = -JhS{S + 1) + JHStotiStot + l)/2 - H.ffSl, - i?/(Mife-2 + f^k+i)/^, (3) 

which determine the total quantum spin number of the spin-S* Heisenberg dimer and its projection 
towards the z-axis, respectively. According to the basic laws of quantum mechanics, the total 
quantum spin number of the spin-S" dimer gains 25" + 1 different values Stot = 0, 1, . . . , 25*, while 
the quantum spin number St^t gains 2Stot + 1 different values S^gf. = —Stot, —Stot + 1, • • ■ , Stot 
for each allowed value of Stot- Then, the energy eigenvalues ([3]) can be straightforwardly used to 
obtain the relevant trace emerging in the last expression on right-hand-side of the equation In 
addition, the resulting expression immediately implies a possibility of performing the generalized 
decoration- iteration transformation [ [371 131] 

Trfc eM-mk) = cMPJhS{S + 1) + pHi{^il^_^ + 

2S n 

X ^ ^ ex.p[—pjHn{n + l)/2] cosh{f3Hejfm) 

n— m— — n 

= Aexp[/3ii/x^fe„2M3fe+i + PHo{^^lk-2 + f4k+i)/2]. (4) 

From the physical point of view, the mapping transformation ^ effectively removes all the inter- 
action parameters associated with a couple of the Heisenberg spins from the fcth diamond unit and 
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replaces them by the effective interaction R and the effective field Hq acting on the remaining Ising 
spins /i3fc-2 and iJ-sk+i only. In this way, one establishes a simple connection between the exact 
solution of the mixed spin-1/2 and spin-S" Ising-Heisenberg diamond chain and the uniform spin- 
1 /2 Ising linear chain with the effective nearest-neighbour interaction R and the effective magnetic 
field Hq. Of course, the transformation relation ([4|) must hold for all possible spin combinations of 
the Ising spins fi3k-2 and fi^k+i- This 'self-consistency' condition then unambiguously determines 
so far not specified mapping parameters A, R and Hq, 

A^cxp[(3JhS{S+1)]{W+W-W^)'^\ /3i? = ln(^^^±^^ , /JHq - - In (^^) , (5) 

which can be uniquely expressed through the functions W± = F{±1) and W = F{0) with 

2S n 

F{x) ^J2Y1 exp[-f3JHn{n + l)/2] cosh[/3m( J/x + Hh)]. (6) 

n— m—~n 

Now, a direct substitution of the transformation ([4]) into the expression ([2]) yields the equality 

Z{f3, Jh, Ji, Hh, Hi) = A^Zo(/3, R, Ho), (7) 

which establishes an exact mapping relationship between the partition function Z of the mixed-spin 
Ising-Heisenberg diamond chain and the partition function Zq of the uniform spin-1/2 Ising linear 
chain with the nearest-neighbour coupling R and the effective magnetic field Hq. It is valuable 
to remark that the mapping relation ([7|) in fact completes our exact calculation of the partition 
function, since the partition function of the uniform spin-1/2 Ising chain can simply be calculated 
within the framework of the transfer-matrix method [ 1461 147] 



Zo(/3,i?,i?o) =exp(iV/3i?/4) 



cosh {(3Ho/2) + y^sinh^ (/3i^o/2) + exp(-/3i?) 



JV 



(8) 



At this stage, the mapping relationship ([7]) between the partition functions can in turn be utilized 
also for a straightforward calculation of the Hclmholtz free energy of the Ising-Heisenberg diamond 
chain 

T = -kBTlnZ = To - TV/cbT In A, (9) 

where !Fq = — fceTlnZo represents the Hclmholtz free energy of the corresponding spin-1/2 Ising 
chain. Subsequently, some other important physical quantities can readily be calculated by using the 
standard thermodynamic relations. Indeed, the sub-lattice magnetization mj and mn reduced per 
one Ising and Heisenberg spin, respectively, can easily be obtained by differentiating equation ([9]) 
with respect to the particular magnetic field acting on Ising and Heisenberg spins 



mn = -^^5i/;^^ = ^^i--o + Soj+^^^+rno+eoj + /:.„^--eoj,(ll) 

In above, the parameters mg and eg denote the single-site magnetization and the correlation 
function between the nearest-neighbour spins of the uniform spin-1/2 Ising chain [ 1461 147'. which 
are unambiguously given by equations ([5]) for the mapping parameters R and Hq through 

1 sinh(/3Ho/2) 
"^0 = 7^ I (12) 
^ y^sinh^ (}3Ho/2) + exp 

cosh mo/2) - Jsinh^ (/3i7o/2) + exp {-pR) 

eo^ml + \ (1 - ml) . (13) 

cosh {PHo/2) + .^sinh^ {(3Ho/2) + exp {-(3R) 
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Finally, the coefficients £^ and Hx emerging in equation (fTTjl mark the expressions £^ = ^ In W± 
and = ^ In whose explicit forms are too cumbersome to write them down here explicitly 
(note that their explicit forms directly follow from the definitions given above after a straightfor- 
ward but rather lengthy calculations). Now, the total magnetization normalized per one spin of 
the diamond chain can be calculated from the expression m = (m/ + 2TO^f)/3. Furthermore, the 
entropy and specific heat of the mixed-spin Ising-Heisenberg diamond chain can directly be calcu- 
lated from the standard thermodynamical relations by computing the first and second temperature 
derivatives of the Helmholtz free energy In this way, one attains the following rigorous results 
for the entropy and specific heat per one site of the investigated Ising-Heisenberg diamond chain 

^ = _LlnZo + ilnA-^[J^5(5+l) + moi?7] 



- - mo + £0 j - ^ (^^ + mo + £0 j - ^ (^2 " j ' 

'^'^ ' - - Too + £o ) + TT^ 1 T + + £0 ) H TT^ ( - - 2£o 

" (15) 



BiVfcB 3 V4 / 3 V4 7 3 V2 





( 9too 


a£o\ 






^ a/3 


dp) 





9£o\ 20^LpdeQ p'^Hidm, 



\ dp dp J 3 dp 3 dp 

where C^^ = ■^lnW± and Cpp — ^^-InW^. As one can see, both afore-listed quantities are 
expressed in terms of the well-known exact results for the partition function Zq, single-site mag- 
netization Too and nearest- neighbour correlation function £o of the corresponding spin-1/2 Ising 
linear chain [ 137] . Our exact calculation of the entropy and specific heat are thus essentially 
completed by evaluating the coefficients £/3, -C^^, £/3/3 and inverse temperature derivatives of 
mo and eo, which are not explicitly given here for brevity. 



3. Results and discussion 



In this part, let us proceed to a discussion of the most interesting results obtained for the 
mixed spin-1/2 and spin-S* Ising-Heisenberg diamond chain. Before doing so, however, let us make 
few remarks on a validity of analytical results presented in the preceding section. It should be at 
first pointed out that all obtained results are rather general as they hold for arbitrary quantum 
spin number S of the Heisenberg spins and also independently of whether ferromagnetic or anti- 
ferromagnetic interactions Jh and Jj are considered. It is also noteworthy that some particular 
cases of the investigated model system have already been examined by the present authors in 
earlier papers [ [Ml [Ml [M] ■ More specifically, the present results reduce to those acquired for 
the Ising-Heisenberg diamond chains with two particular spin values 5 = 1/2 and 1, which have 
been undertaken a rather detailed analysis in references [ [32l [33l [34] ■ Furthermore, the present 
model also contains the mixed spin-1/2 and spin-3/2 Ising-Heisenberg diamond chain as another 
particular limiting case, which has been a subject matter of our preliminary report [[35] revealing 
a series of intermediate plateaux in the magnetization process of this quantum spin chain. 

With this background, we will restrict ourselves in this paper to an analysis of another particular 
spin cases of the mixed-spin Ising-Heisenberg diamond chain with the aim to shed light on how the 
magnetic behaviour of the model under investigation depends on the quantum spin number S of 
the Heisenberg spins. Namely, this systematic study should provide a deeper understanding of the 
role of quantum fluctuations in determining the overall magnetic behaviour, because the quantum 
effects should become less significant by increasing the quantum spin number S. As we are mainly 
interested in the examination of the geometric frustration, we will assume in what follows the 
antiferromagnetic character of both interaction parameters Jh > and Jj > 0. Next, it is also 
convenient to rescale all interaction parameters with respect to the Ising- type interaction J/, which 
will henceforth serve as the energy unit. In doing so, one effectively introduces the following set 
of dimensionless parameters: a = Jh/Ji, h ~ Hi/ Ji ~ Hh/Ji, and t = k^T/Jj, as describing a 
strength of the geometric frustration, the external magnetic field, and the temperature, respectively. 
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3.1. Ground-state properties 

Let us begin our discussion by considering possible spin arrangements to emerge in the ground 
state of the mixed spin-1/2 and spin-S Ising-Heisenberg diamond chain. For this purpose, the 
ground-state phase diagrams in the a — h/S plane are displayed in figure [2] for four particular spin 
cases. As one can see from figures [2][a)-(c), the mixed-spin diamond chain with finite quantum spin 
numbers S — 3/2, 2, and 5/2 have quite similar ground-state phase diagrams, which differ mainly 
in the total number of possible ground states. Apart from the semi-classically ordered ferrimagnetic 
phase (FRI), saturated paramagnetic phase (SPP) and frustrated phase (FRU), there also appear 
25* — 1 quantum ferrimagnetic (to be denoted as QFIi, QFI2, . . . , QFl2s_i) and 25—1 quantum 
ferromagnetic (to be denoted as QFOi, QFO2, . . . , QFO2S-1) phases. Spin arrangements of the 
relevant phases can be unambiguously characterized by means of the following eigenfunctions and 
single-site magnetization 

N N 

|FRI) = l[\-hk-2l[\S,S)3k^i^sk, 

k=l k=l 

mi = -1/2, niH = S, m/msat = 1 - 2/(45 + 1); (16) 

N N 

3fc-l, 3k, 

k=l k=l 

mi = 1/2, mn = 5, m/msat = 1; (17) 

N N S 

fe=l fe=l + ^ l=-S 

mi = 0, mn = 0, m/mgat = 0; (18) 

N N S 

fe=l k=l l=S-j 

mi = -1/2, mH = S- j/2, m/msat = 1 - 2(j + l)/(45 + 1); (19) 

N N S 

IQFO^.) = ni+)3^-2n E {-lf^'Ai^2S-l~,\l,2S~l-j)3k-l,3k, 
k=l k=l l=S-j 

mi = 1/2, mH = S- j/2, m/msat = 1 - 2^/(45 + 1); (20) 

where j = 1, 2, . . . , 25—1. Note that the first product in the afore- listed eigenfunctions is taken over 
all Ising spins (|±) denotes /i^ = ±1/2), while the second one runs over all Heisenberg dimers. The 
coefficients Ai^2S-i-j emerging in the last two eigenfunctions and (|20p represent probability 
amplitudes for finding the Heisenberg spin pairs in the spin state \l,2S — I — j) and these are listed 
for several particular spin cases in the Appendix together with the complete analytic form of the 
eigenfunctions QFIj and QFO^ . Finally, m^at labels the saturation magnetization normalized per 
one spin of the diamond chain. 

It is quite evident from the set of equations p6)) " ([20l) that first two phases FRI and SPP exhibit 
spin arrangements, which are commonly observed also in the pure Ising systems. By contrast, the 
quantum entanglement between spin states of the Heisenberg spin pairs represents an inherent 
feature of the phases FRU, QFIj and QFOj, which cannot be observed in the pure Ising systems. 
According to their location in the ground-state phase diagrams, the phases with entangled spin 
states of the Heisenberg spins appear as a result of a competition between the interaction pa- 
rameters a and h, i.e., owing to a mutual interplay between the geometric frustration generated 
by the competition between the Heisenberg- and Ising- type interactions Jh and J/, respectively, 
the quantum fluctuations arising from the Heisenberg interaction Jh, and the effect of applied 
magnetic field H. It is also noteworthy that the phases QFI^ and QFOj basically differ from each 
other just in a relevant spin alignment of the Ising spins, which are oriented antiparallel (parallel) 
with respect to the total spin of the Heisenberg dimers in QFIj (QFOj). Hence, the overall spin 
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arrangement of the phases QFlj (QFOj) has typical features of the quantum ferrimagnetic (ferro- 
magnetic) phase with a significant quantum reduction of the sublattice magnetization rriH closely 
related to a quantum entanglement of the Heisenberg spins. 

A profoundly different situation emerges when one is assuming the limit of classical Heisenberg 
spin S —^ (X), whose ground-state phase diagram is shown in figure[2l^d). As one can see, solid lines 
separate just four different phases, namely, the FRI phase, the SPP phase and two phases denoted 
as {FRIi, FRI2, FRIoo} and {FROi, FRO2, FROoo}- As could be expected, the spin 
arrangements emerging in the phases FRI and SPP are very analogous to the ones of the diamond 
chains with finite values of the Heisenberg spins. In the FRI phase, one actually finds rn/ = —1/2, 




mn — > 00, m/nisat — 1, which indicate the semi-classical ferrimagnetic character of this phase. 
In the SPP phase being stable at sufficiently strong fields, one observes mj ~ 1/2, tuh 00, 
m/msat — 1 implying that all Ising as well as Heisenberg spins are oriented in a direction of the 
applied external magnetic field. A detailed examination of another two possible phases {FRIi, FRI2, 
. . . , FRIoo} and {FROi, FRO2, . . . , FROoo} reveals that the former (latter) phase has ferrimagnetic 
(ferromagnetic) character, since the sublattice magnetization m/ = —1/2 (1/2) is oriented in the 
opposite (same) direction as the sublattice magnetization mn, which linearly increases with the 
applied magnetic field h on account of a gradual rotation of the classical Heisenberg spins into the 
external-field direction. The linear increase of the sublattice magnetization mn can alternatively be 
viewed also as a smooth sequence of infinite number of phase transitions between the phases FRIi , 
FRI2, . . ., FRIoo (or FROi, FRO2, . . ., FROoo) representing classical analogs of the quantum phases 
QFIi, QFI2, . . ., QFIoo (or QFOi, QFO2, . . ., QFOoo)- However, each phase transition between 
those phases is connected with just an infinitesimal increase of the sublattice magnetization ttih if 
S ~* 00 and hence, these field-induced transitions cannot be regarded as true phase transitions. 
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3.2. Magnetization process 

Now, let us turn our attention to the magnetization process of the investigated model system. 
It is useful to mention, however, that possible magnetization scenarios for the Ising-Heisenberg 
diamond chain with three lowest possible spin values of the Heisenberg spins S = 1/2, 1 and 3/2 
have already been particularly examined in our previous papers [ I34[ I35j to which the interested 
reader is referred to for more details. In this respect, we merely depict in Fig. [3] the total and 
sublattice magnetization for another particular case S = 5/2, which should bring insight into how 
possible magnetization scenarios evolve by increasing significantly a magnitude of the Heisenberg 
spins. One may immediately come to the following general conclusions valid for the Ising-Heisenberg 
diamond chains with the arbitrary finite value S of the Heisenberg spins by combining the formerly 
published results [see figures 2(a), 3, 7(a), 8 in reference [ [34| and figure 1 in reference [[35]] 
with the results shown in figures [2lja)-(c) and [3l there always exists the critical value of the 
frustration parameter ac = (25') ~^ below which just a single intermediate magnetization plateau 
appears in the magnetization curve and above which more diverse magnetization scenarios can 
be in principle observed with two, three,. . 2S intermediate magnetization plateaux. The zero- 
temperature magnetization curves corresponding to the particular case with a < ac essentially 
reflect the field-induced phase transition FRI-SPP, while they might reflect another five different 
sequences of phase transitions if a > etc, namely, 

1. QFIj-QFIj_i- • • • -QFIi-FRI-SPP {j = 1,2,. . . ,25-1) 

2. FRU-QFI25-2-QFI25-3 QFIi-FRI-SPP 

3. FRU-QFO2S-1-QFO2S-2 QFOj-QFIj_2-QFIj_3 QFIi-FRI-SPP (j = 3,4,. . . ,25-1) 

4. FRU-QFO2S-1-QFO2S-2 QFO2-FRI-SPP 

5. FRU-QFO2S-1-QFO2S-2 QFOi-SPP 

Moreover, it is also worthwhile to remark that all intermediate plateaux identified in the relevant 
field dependences of the total magnetization appear at 1/(45'+ 1), 3/ (45'+ 1), . . . , (45 — 1)/(45'+ 1) 
of the saturation magnetization. The closer mathematical analysis reveals that all these fractional 
values satisfy the Oshikawa-Yamanaka- Affeck rule p {Su — m) £ N [ 48 , which has been proposed 
as a necessary condition for the formation of quantized plateaux {p is a period of the ground state, 
Su and TO are the total spin and total magnetization of the elementary unit). 

To compare the magnetization scenario of the diamond chain with the finite and infinite value 
of the Heisenberg spin, typical magnetization versus magnetic field dependences are depicted in 
figure m for two particular spin cases 5 = 4 and S ^ 00. Note that the magnetization curves 
corresponding to the limit of classical Heisenberg spin S —^ 00 are for clarity shown as broken 
lines. It is easy to observe from this figure that the diamond chain with 5 — s- 00 has qualitatively 
different magnetization process in comparison with its quantum counterpart with the finite (albeit 
very high) value of the Heisenberg spins. Actually, the magnetization curves display a linear increase 
of the magnetization with increasing magnetic field until the magnetization reaches its saturation 
value if one considers the classical limit 5' — > 00 instead of the sharp stepwise magnetization 
curves, which appear for any finite value of the Heisenberg spin. The observed linear increase 
of the magnetization in the phases {FRIi, FRI2, ...,FRIoo} and {FROi, FRO2, FROoo} 
is consistent with the fact that the intermediate magnetization plateaux gradually shrink as the 
quantum spin number 5 increases until they entirely merge together into the linear magnetization 
vs. external field dependence in the classical limit 5 — > cxd. 

3.3. Specific heat 

In this part, let us take a closer look at the specific heat versus external field dependences. 
Some typical variations of the specific heat as a function of the dimensionless magnetic field di- 
vided by the Heisenberg spin 5' are depicted in figure[5]for the particular spin case 5' ~ 5/2, several 
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Figure 3. The total and sublattice magnetization against the external magnetic field at zero 
temperature t = 0.0 for the particular spin case 5* = 5/2 and the fi-ustration parameter: (a) 
a = 0.0, (b) a = 0.75, (c) a = 1.25, (d) a = 1.6, (e) a = 1.78, and (f) a = 2.0. 



values of the frustration parameter a and three different temperatures. To enable a direct compar- 
ison, the frustration parameter a is chosen so as to match four different magnetization scenarios 
plotted in figures [3] (b) , (d)-(f). The displayed sets of C{h) curves thus basically reflect four mag- 
netization processes, which are accompanied by the following sequences of the field-induced phase 
transitions: QFI4-QFI3-QFI2-QFI1-FRI-SPP [figure El^a)], FRU-QFO4-QFI2-QFI1-FRI-SPP [fig- 
ure [^b)], FRU-QFO4-QFO3-QFO2-FRI-SPP [figure Eljc)], and FRU-QFO4-QFO3-QFO2-QFO1- 
SPP [figure [Sfd)] . It can be clearly seen from these figures that the thermal trend of the specific 
heat as a function of the magnetic field is quite similar for all investigated magnetization scenarios: 
C{h) curves have rather irregular shapes with several broad maxima at higher temperatures (see 
the curves labeled as i = 0.6) that develop into more pronounced peaks gradually moving towards 
transition fields as the temperature is lowered (see for instance the curves t = 0.2). Naturally, the 
total number of peaks in the low-temperature C{h) curves depends on the value of the frustration 
parameter a. More specifically, the zero-temperature magnetization process formed by the sequence 
of field-induced transitions between the quantum ferrimagnetic phases [such as QFI4-QFI3-QFI2- 
QFIi-FRI-SPP shown in figure El^a)] characterize at sufficiently low temperatures identical double 
peaks in C{h) curves being symmetrically centered around the respective transition fields. On the 
other hand, one additional single peak is always detected in the low-temperature C{h) curves in 
the limit of vanishing external field whenever the frustration parameter a drives the system into 
the disordered FRU ground state [see the case t = 0.05 in figures[5l^b)-(d)]. The origin of this single 
peak lies in thermal excitations of the frustrated Ising spins, which have tendency to align towards 
the external-field direction at low enough temperatures. Another interesting observation coming 
from figure is that almost all peaks are identical both in their height as well as width at low 
temperatures. The only exceptions constitute double peaks located around the critical fields asso- 
ciated with the phase transitions QFO4-QFI2 and QFO2-FRI [see insets in figures [S^b) and (c)], 
whose origin significantly differs in some respects from the origin of others. As a matter of fact, 
both afore-mentioned double peaks result from intensive thermal excitations of Ising as well as 
Heisenberg spins, because both of them are changing their spin state at the transition fields be- 
tween the phases QF0j-QFIj_2 (j = 1, 2, . . . , 25* — 1) and QFO2-FRI. Contrary to this, the other 
peaks reflect vigorous thermal excitations of only the Ising spins or only the Heisenberg spins as it 
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Figure 4. The total magnetization reduced with respect to its saturation value versus the external 
field at zero temperature t = 0.0 for the spin cases 3 = 4 and 5* cxd and the frustration 
parameter: (a) ct = 0.05, (b) a = 0.75, (c) a = 1.7, and (d) a = 2.0. 



directly follows from a detailed investigation of the low-temperature magnetization process shown 
in figure [31 

At this point, it is worthwhile to compare our results with those obtained by Efremov and 
Klemm for the antiferromagnctic spin-^* Heisenberg dimer [ 1491 150] . These authors have found 
that the low-temperature C{h) curves of the antiferromagnetic spin-^ Heisenberg dimer exhibit a 
quite universal dependence with several marked double peaks symmetrically centered around level- 
crossing fields, which arise from thermal excitations between the ground state and first excited state. 
The relevant expression for the specific heat substantially simplifies in a close vicinity of the level- 
crossing fields at low temperatures, where the specific heat is well approximated by the equation 
(36) of the reference [[50], because the only significant contribution to the specific heat comes 
from thermal excitations between the ground state and the first excited state. In this respect, the 
specific heat goes exponentially to zero as one reaches the level-crossing field and there also appear 
two symmetrically centered peaks around the level-crossing field whose heights are given by the 
condition C^'^'^/Nk^ = (c/coshc)^ « 0.439229 (c « 1.199679 is the solution of the transcendent 
equation tanh c = 1 / c) . It is noteworthy that the specific heat of the mixed-spin Ising-Heisenberg 
diamond chain is obviously governed by the same asymptotic expression near the transition fields 
as evidenced by the height and position of the relevant peaks in the field-dependence of the specific 
heat. In addition, it is also quite interesting to mention that the single peak, which is detected in 
the low-field tail {h — > 0) of the specific heat on assumption that the Ising-Heisenberg diamond 
chain is driven by the sufficiently strong frustration a > 1 towards the disordered FRU phase [see 
figures[5Kb)-(d)], is of the same height as the symmetrically centered double peaks. It is worthwhile 
to recall that this single peak originates from the field-induced splitting of the highly-degenerate 
lowest energy level and the same height of this peak indicates considerable contribution of the single 
spin-flip excitations of the Ising spins to the specific heat. On the other hand, the nonuniform double 
peaks observed in low-temperature C{h) curves of the mixed-spin Ising-Heisenberg diamond chains 
cannot be described by the simple formula proposed by the analysis of the antiferromagnetic spin-S' 
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Figure 5. The specific heat versus the external magnetic field at various temperatures t — 0.05, 
0.2, 0.6 for the spin case S = 5/2 and the frustration parameter: (a) a = 0.75, (b) a = 1.6, (c) 
a = 1.78, and (d) a = 2.0. The details show the non-symmetric double peaks in the enlargened 
scale (for details see the text). 



Heisenberg dimer ["SD]. This implies that thermal excitations between more than two energy levels 
might possibly come into play at these particular transition fields, where both the Ising spins as 
well as the Heisenberg spins change their spin state. 



3.4. Enhanced magnetocaloric effect 

Finally, we will turn our attention to an investigation of the adiabatic demagnetization studied 
in connection with the enhanced magnetocaloric effect. For this purpose, the most interesting 
results for the field variations of the temperature are plotted in figure [6] for four different values of 
the Heisenberg spins S by keeping the entropy constant at the particular value of the frustration 
parameter a = 2.0. Note that this particular value of the frustration parameter was chosen so that 
the highly-degenerate FRU phase will constitute the ground state in the limit of vanishing external 
magnetic field. Namely, it has been recently proved that one achieves a substantial enhancement 
of the cooling rate during the adiabatic demagnetization of highly frustrated spin systems in 
comparison with the adiabatic demagnetization of unfrustrated spin systems [ l5Tl[52l[53l[54l[55ll56j . 
As one can readily see from figures [5l[a)-(d), the mixed-spin Ising-Heisenberg diamond chains 
exhibit a pronounced valley-peak structure in the field dependence of the temperature at a fixed 
value of the entropy, which differ mainly in the total number of peaks (valleys) that gradually 
increases with the Heisenberg spin S. The most obvious drop (grow) of the temperature can 
always be found in the vicinity of zero field and transition fields at which the system undergoes 
zero-temperature phase transitions (see figure [2] for the ground-state phase diagram) . It can be 
also observed from figure [6] that a cooling down to the lowest temperatures is achieved only if the 
entropy is less than or equal to the residual entropy ares = Sres/'^N = ^ln2 of the disordered 
FRU phase, otherwise, the nonzero temperatures are finally reached as the external magnetic field 
vanishes. 
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Figure 6. Adiabatic demagnetization in the form of temperature versus external magnetic field 
dependence for several values of the entropy per one spin a = S/3N, the frustration ratio 
a = 2.0 and the Heisenberg spins: (a) 5" = 3/2, (b) S = 2, {c) S = 5/2, and (d) S = 4. For 
better orientation, broken curves depict the relevant dependences when the entropy is equal to 
the residual value ares = | ln2 = 0.231049. 



To analyze the cooling rate of the adiabatic demagnetization in a more detail, the adiabatic 
magnetocaloric rate (dt/dh)^ is depicted in figure [7] against the external magnetic field for the 
particular spin case 5 = 5/2 and the frustration ratio a = 2.0 by keeping the entropy per one spin 
(T = S/3N constant at four different values. As one can see, the most rapid cooling and heating 
rate can be observed during the adiabatic demagnetization when the entropy is close enough to 
the residual entropy ares = |ln2 of the disordered FRU phase. In the particular case a = ares, 
one even observes asymptotically infinitely fast cooling (heating) as some particular transition field 
is approached from above (below). It is worthy of notice here that a relatively fast cooling of the 
frustrated diamond chain, which is observable in a certain range of fields and temperatures, might 
be of practical importance for a cooling purpose. The cooling effect is of technological relevance 
only if the relevant cooling rate exceeds the one of paramagnetic salts {dt / dh)^^'^'^ — t/h and the 
applied external magnetic fields are within the experimentally accessible range. From this point 
of view, the adiabatic demagnetization of the frustrated diamond chains becomes fairly efficient 
in producing the enhanced magnetocaloric effect just if the entropy is chosen close enough to 
its residual value ares and the applied magnetic field is sufficiently small. This makes from the 
frustrated diamond chain systems promising candidates for being efficient refrigerant materials, 
which would allow reaching ultra-low temperatures unlike the paramagnetic salts that usually 
exhibit a spin-glass transition. However, the cooling rate basically diminishes when the entropy 
is selected far from the residual value ares and this may represent rather inconvenient property 
hindering from real-world applications by reaching the ultra-low temperatures with the help of the 
adiabatic demagnetization. 
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Figure 7. Adiabatic magnotocaloric rate {dt/dh)„ as a function of the external magnetic field 
for the particular case with the Heisonbcrg spins S = 5/2, the frustration ratio a = 2.0, and 
four different values of the entropy per one spin (cr = S/3N): (a) a = 0.02, (b) a = 0.10, (c) 
ares = I ln2 = 0.231049, and (d) a = 0.30. 



4. Conclusions 



In the present article, the mixed spin-1/2 and spin-5 Ising-Heisenberg diamond chain has ex- 
actly been solved by combining three precise analytical techniques: the Kambc projection method, 
the generalized decoration-iteration transformation, and the transfer-matrix method. Within the 
framework of this exact analytical approach, the ground-state phase diagrams and the magnetiza- 
tion process were particularly examined in dependence on a strength of the geometric frustration, 
as well as, the quantum spin number S of the Heisenberg spins. In addition, the specific heat as a 
function of external magnetic field has been explored in detail along with the adiabatic demagne- 
tization process, which was examined in connection with the enhanced magneto caloric effect. 

The main goal of the present work was to shed light on how the magnetic behaviour of the 
mixed-spin diamond chains depends on a magnitude of the Heisenberg spins. The rather systematic 
study of the Ising-Heisenberg diamond chain models with different values of the Heisenberg spins 
indeed brought a deeper insight into how a diversity of the magnetization scenarios is related 
to the magnitude of the Heisenberg spins. It has been shown that the mixed spin-1/2 and spin- 
s' Ising-Heisenberg diamond chains generally exhibit multistep magnetization curves with up to 
25 intermediate magnetization plateaux provided that there is a sufficiently strong geometric 
frustration. Beside this, we have also performed the comparative study of the classical limit of 
the Ising-Heisenberg diamond chain with the classical vector spins S ^ oo, which enabled us to 
discern the typical quantum features (e.g. multiply steps in the magnetization curves) from the 
typical classical ones (e.g. the linear magnetization vs. external field dependence). 
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Appendix 

This appendix explicitly enumerates the coefficients Ai^2S-i~j: which emerge in the eigenfunc- 
tions (|19p and (pn|) characterizing the quantum phases QFlj and QFOj (j = 1, 2, . . . , 25* — 1), that 
determine probability amplitudes for finding the Heisenberg spin pairs in the spin state \l, 2S—l—j). 
In what follows, the relevant probability amplitudes are listed for three particular spin values S* = 1, 
3/2 and 2 of the Heisenberg spins together with the explicit form of the phases QFIj and QFOj. 
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